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ABSTRACT 



We investigate the accretion process in high-luminosity AGNs (HLAGNs) 
in the scenario of the disk evaporation model. Based on this model, the thin 
disk can extend down to the innermost stable circular orbit (ISCO) at accretion 
rates higher than 0.02MEdd! while the corona is weak since part of the coronal 
gas is cooled by strong inverse Compton scattering of the disk photons. This 
implies that the corona cannot produce as strong X-ray radiation as observed 
in HLAGNs with large Eddington ratio. In addition to the viscous heating, 
other heating to the corona is necessary to interpret HLAGN. In this paper, we 
assume that a part of accretion energy released in the disk is transported into 
the corona, heating up the electrons and thereby radiated away. We for the first 
time, compute the corona structure with additional heating, taking fully into 
account the mass supply to the corona and find that the corona could indeed 
survive at higher accretion rates and its radiation power increases. The spectra 
composed of bremsstrahlung and Compton radiation are also calculated. Our 
calculations show that the Compton dominated spectrum becomes harder with 
the increase of energy fraction (/) liberating in the corona, and the photon index 
for hard X-ray (2 — 10 keV) is 2.2 < T < 2.7. We discuss possible heating 
mechanisms for the corona. Combining the energy fraction transported to the 
corona with the accretion rate by magnetic heating, we find that the hard X-ray 
spectrum becomes steeper at larger accretion rate and the bolometric correction 
factor (Lboi/-^2-iokcv) increases with increasing accretion rate for / < Jr, which 
is roughly consistent with the observational results. 

Subject headings: accretion: accretion disk-galaxies: active-X-rays: galaxies 
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1. INTRODUCTION 

Accretion of gas onto the central supermassive black hole is one of the fundamental 
astrophysical processes responsible for high-efficient energy release in active galactic nuclei 
(AGNs). According to the difference in the luminosity, AGNs are classified into two types, 
low-luminosity AGNs (LLAGNs) and high-luminosity AGNs (HLAGNs). The LLAGNs 
show hard power-law X-rays and/or radio jets; while the HLAGNs are characterized by 
the big blue bump around UV waveband, soft X-ray excess, Fe K a lines at about 6.4 keV, 
the power-law X-rays with photon index V ~ 1.9 and cutoff at about a few of 100 keV. 
It is speculated that the different characteristics of LLAGNs and HLAGNs are caused 
by different accretion mechanism at different accretion rates, in a similar way to black 
hole X-ray binaries. The accretion in LLAGNs is via an advection dominated accretion 
flow (ADAF) (Ichimaru 1977; Narayan & Yi 1994, 1995a, 1995b), in which most of the 
accretion energy is stored in the gas as entropy and only a small fraction is radiated in hard 
X-rays(see Kato, Fukue & Mineshige 2008 for a review). While in HLAGNs, the common 
accepted accretion model is a hybrid accretion disk, i.e., a cool disk is embedded in a hot 
corona. The optical to UV emission is from the optically thick and geometrically thin 
accretion disk (e.g., Shakura & Sunyaev 1973; Lynden-Bell & Pringle 1974) and the hard 
X-ray is from Compton scattering of the soft photons from the disk by the hot electrons in 
the optically thin corona with T c ~ 10 9 K. 

Observations show that the X-ray luminosity contributes a large fraction to the 
bolometric luminosity in HLAGNs. This can be seen from individual spectra (e.g. Elvis 
et al.1994; Ho 2008). The hard X-ray bolometric correction, which is defined as the 
ratio of bolometric luminosity to 2 — 10 keV luminosity, is typically around 40 — 70 with 
Eddington ratio above 0.1, while below 0.1 it is typically 15 — 25 (Vasudevan & Fabian 
2009), which confirms the importance of X-ray radiation in LLAGNs. Fits to the X-ray 
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spectra of HLAGNs by Comptonization model yield an electron temperature of ~ 10 9 K in 
the optically thin hot corona (e.g. Liu et al. 2003). A question arises as, how can the gas 
in the corona retain at so high a temperature when continuously radiates strong X-rays in 
HLAGNs? Where is the heating energy from? 

Theoretically, gas accreted from interstellar medium or secondary is cold. In the way 
of accretion to the central black hole, some of the gas evaporates to the corona and then 
accretes inward in the corona. In HLAGNs, gas is dominantly accreted via a thin disk, 
only a small fraction of gas is accreted in the hot corona. Therefore, the energy released 
in the corona through accretion is much smaller than that in the disk. This implies that 
the radiation from the hot corona is much weaker than that from the disk in a steady 
disk-corona flow. With strong radiation in X-rays as observed in HLAGNs, the corona is 
over-cooled in a thermal timescale and collapses to the disk unless there is other energy 
supply to the corona besides the viscous heating. 

Recent investigations on the spectral energy distribution (SED) in HLAGNs found 
that the bolometric correction increases with increasing Eddington ratio Lboi/^Edd 
(-^Edd — 1-26 x lO 46 (M B H/lO 8 M )ergs s" 1 ), whereas independent on the black hole mass 
(Vasudevan & Fabian 2007, 2009; Lu & Yu 1999; Wang et al. 2004; Liu et al. 2009b). 
This correlation provides clues and constraints to possible heating mechanisms to the disk 
corona in HLAGNs. 

In this work, we study the corona flow at high accretion rates in the frame of disk 
corona evaporation/condensation model. By inclusion of additional heating to the corona, 
we investigate the detailed interaction, i.e. the mass and energy exchange between disk and 
corona, we self-consistently determine the corona structure and disk radiation features. The 
spectrum from the disk and corona is calculated by Monte Carlo simulations. Comparing 
the spectrum with observational features, we discuss the possible heating mechanism in 
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HLAGNs. Our aim is to interpret the SED and its correlation with Eddington ratio. 

HLAGNs include radio-loud and radio-quiet Seyfert Is and QSOs. Since the radio-loud 
AGNs are assumed to be powered by the jet and the beaming effect will dilute the X-ray 
radiation from the central accretion engine, here we focus on accretion mechanism for the 
radio-quiet AGNs. Most narrow line Seyfert Is are also radio-quiet HLAGNs. However, 
since they accrete at rate near Eddington accretion rate and the character of their spectra 
(e.g. soft X-ray excess and steeper/softer hard X-ray emission) are dramatically different 
from other HLAGNs, they are not include here and will be investigated in the future work. 

The paper is organized as follows. In section 2 we present the model. The numerical 
results are shown in the section 3. Section 4 and Section 5 are the discussion and conclusion, 
respectively. 

2. MODEL 

2.1. Conception of the Model 

The disk evaporation model was proposed by Meyer & Meyer-Hofmeister (1994) for 
dwarf novae and developed to black hole binaries by Meyer, Liu & Meyer-Hofmeister 
(2000a). In this model, an optically thin, hot corona is assumed to lie above a thin disk. 
Vertical conduction is important since the disk temperature is much lower than that of the 
corona. As the corona radiates inefficiently, the conductive heat accumulates downwards 
and at the transition layer it heats up some cool gas until an equilibrium density is reached 
so as to radiate completely the conductive flux. The evaporated gas is then accreted 
through the corona to the central black hole. Such a corona is heated up by viscosity and 
the accretion is supplied by mass evaporation from the disk. Further investigations on 
effects of decoupling between ions and electrons and Compton cooling in the corona (Liu 
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et al. 2002a), the magnetic pressure (Meyer- Hofmeister & Meyer 2001; Qian, Liu & Wu 
2007), the change of viscosity parameter (Qiao & Liu 2009) and the possible condensation 
are included. The model has been successfully used to interpret the state transition and 
observational features in low/hard state of the black hole binaries. As the disk corona 
is independent on the black hole mass, it is also applied to AGNs in explaining different 
types of spectra and the absence of the broad line region in LLAGNs (Liu et al. 2009a). 
Nevertheless, when it is compared with observations in HLAGN, we find that in the corona 
strong Compton cooling leads to the corona so weak that it cannot produce the hard X-ray 
as high as observed. There must be other energy sources besides the viscous heating for the 
corona. 

There are some pioneer works assuming that a certain fraction of local gravitational 
energy is liberated directly in the corona (e.g. Nakamura & Osaki 1993; Haardt & Maraschi 
1991, 1993), particularly in the frame of a magnetized disk corona (e.g., Di Matteo 1998; 
Di Matteo, Celotti, & Fabian 1999; Miller & Stone 2000; Merloni & Fabian 2002; Liu et 
al. 2002b, 2003; Kawanaka, Kato, & Mineshige 2008; Cao 2009). In these investigations 
it is demonstrated that the corona can indeed produce as strong X-ray as observed in 
HLAGNs if a large fraction of accretion energy is released in the corona. Nevertheless, the 
mass exchange between the two-phase accretion flows was not taken into account. Here we 
consider both radiation and mass coupling between disk and corona and study the detailed 
vertical structure of the coronal flow above a thin disk supposing some fraction (/) of the 
gravitational energy is added into the corona. 
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2.2. Physical Processes and Differential Equations 

The corona can be described by the following equations. Equation of state 

^ = ^(Ti + T c ), (1) 

where \x = 0.62 is the molecular weight assuming a standard chemical composition 
(X = 0.75, Y = 0.25) for the corona. Here we assume n\ = n c , which is strictly true only 
for a pure hydrogen plasma. 



Equation of continuity 



where = — is the radial component of velocity, V s = (-^) 1 ^ 2 is the isothermal sound 
speed, v v = ^-(1 + ^-)~ 4//3 is the angular component of velocity. Here, the partial 
derivative of mass flux with respect to radius is approximated as, 9 R J^ R ~ —Vm^PVr, 
where rju depends on the net mass gain/loss rate through the radial boundaries (Meyer- 
Hofmeister & Meyer 2003; Liu et al. 2004). In this work, rju = 1 is taken for the case 
that no mass enters through the outer boundary and all the mass flowing in the corona is 
contributed by the evaporation (Meyer et al. 2000a). 

Equation of the ^-component of momentum is 

dv z dP GMz 

In the hot corona, because of heavier mass of ions than that of electrons, the viscous 
heating raises the ion temperature first and ions are cooling by coulomb collision with 
electrons and both radial and vertical advection. Here we don't include direct heat to 
electrons. The ion thermal conduction is not taken into account because of the long mean 
free path compared to the vertical scale height of the corona. So the energy equation for 
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where ?7e is the energy modification parameter, ?7e = %i + 0.5. The difference 
between 7/e and 7/m comes from the derivative of energy flux with respect to radius, 
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its derivative « — -| since the potential, kinetic, and thermal specific energies all scale 
approximately as 1/R (Meyer-Hofmeister & Meyer 2003; Liu et al. 2004). In equation (jlj), 
qi e is the exchange rate of energy between electrons and ions through coulomb collision and 
is described as 

/9\5Sm. UT.5 

(5) 
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where m ; and m e are the proton and electron masses, k is the Boltzmann constant, c is 
the light speed, or is the Thomson scattering cross section and In A = 20 is the Coulomb 
logarithm. 

Different from the previous works, besides the viscous heating, a fraction of local 
gravitational energy is assumed to directly heat the electrons in corona, i.e., 

3GMM 
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add 



8ttR 3 



[1 - (SRs/R) 1 ' 2 } x /, 



(7) 



where M is the total accretion rate in the disk corona. The physical meaning of this 
additional heating term will be discussed in the section of Discussion. To simplify the 
calculations, this additional heating flux is supposed to distribute along the vertical 
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direction in the form similar to that of Compton cooling, i.e. 
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(9) 



and A is determined by the integration along z from the lower boundary to upper boundary, 
J* 1 Qa.dddz = Q a dd, which gives the relation of / and A, 
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Therefore, the energy equation for both the ions and electrons is 
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In this equation, n e niL(T e ) is the bremsstrahlung cooling rate and F c is the thermal 
conduction (Spitzer 1962), 

5 dT 

F c = -K Ti-^, (12) 

az 

with kq = 10~ 6 ergs~ 1 cm~ 1 K~5 for fully ionized plasma. 



For compton cooling rate, gcmp, we mainly consider inverse Compton scattering of the 
soft photons from the disk, while the ones contributed from the bremsstrahlung cooling is 
not included. 

eff (13) 



AkT, 

<?Cm P = ^n e a T c 



m e c* 2 

where a the radiation constant. T e g is the effective temperature of the underlying thin disk, 
fulfilling 

(14) 



87Ti? 3 
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where Md is the mass accretion rate in the thin disk, which depends on the distance because 
of evaporation, 

M d (R)=M-M evap (R), (15) 

with M evap (i?) the integrated evaporation rate from the outer edge of the disk to the distance 
R, which can be approximated by the one- zone evaporation rate, M evap = 2nR 2 (pv z ) ( see 
Liu et al. 2002a). The negative term associated with —fM in Eq. fTHl) corresponds to the 
transportation of accretion energy from the disk to the corona, which results in a lower 
temperature in the disk. Here, we have not considered back reaction (heating of a thin disk 
by coronal illumination) and its justification will be discussed later. 

These five differential equations, Eqs. ([2]), ([3]), (@J, ffTTl) . f fT2l . which contain five 
variables P(z), Ti(z), T e (z), F c (z), and pv z , can be solved with five boundary conditions. 



2.3. Boundary Conditions and Numerical Method of Computation 

At the lower boundary, the temperature of the gas should be the effective temperature 
of the accretion disk. Liu et al. (1995) showed that the coronal temperature increases from 
effective temperature to 10 6 ' 5 K in a very thin layer of nearly constant pressure. Its physics 
is described by the balance between thermal conduction and radiation loss. So a relation 
can be established between temperature and heat flux which can be scaled according to the 
pressure. Combining with Shmeleva-Syrovatskii relation as F c = — 1 \ 1 p (Shmeleva 

& Syrovatskii 1973), the lower boundary conditions can be approximated (Meyer et al. 
2000a) as, 

T i = T c = 10 6 - 5 K, 

(16) 

F c = -2.73 x 10 6 P at z = z , 
There is no pressure and no heat flux at infinity, which requires sound transition at 
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some height z — z\. So we constrain the upper boundary as, 

F c = and v z = V s at z = Z\. (17) 

With such boundary conditions, assuming an initial value A and a pair of lower boundary 
values for P and (pv z )o, we start the integration along z. If the trial values for P Q and 
(pv z )o fulfill the upper boundary conditions, the initial value of P$ and (pv z )o can be taken 
as the solutions of the differential equations. The value of / can then be calculated out from 
Eq.f llOp . The results for a series of / = 0.0, 0.1, 0.2, 0.3, 0.4 are given in the following 
section. When / = 0.0, which means that there are no other energy heating for the corona 
except the viscous heating, this result is similar to the former works and will be compared 
with the results for / > 0.0. 



2.4. Spectra from the Disk and Corona 

With the detailed calculations of the corona structure, the temperature and density in 
the corona and the effective temperature are determined. We then calculate the spectrum 
from such a disk and corona for given structure. 

In the corona, the emissions are from bremsstrahlung radiation and inverse Compton 
scattering of disk photons. For bremsstrahlung cooling spectrum, we adopt the Eddington 
approximation and two-stream approximation (Rybicki & Lightman 1979; Manmoto et al. 
1997) to calculate the radiative flux, 

F v = —B v [l - exp(-2v / 3r;)], (18) 

here B v is Planck function, r* = ^k v H is the vertical absorption optical depth. Assuming 
that the flows is local thermal equilibrium, then k v = Xv/4ffB v , wherex^ = X^.brcms is 
emissivity of bremsstrahlung (Narayan & Yi 1995). 
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Since we focus on HLAGNs, of which the accretion rate is very high and the disk is 
the dominant accretion flow, the soft photons for inverse Compton scatter are supplied 
by the thin disk blackbody radiation. However, the ones from bremsstrahlung emission 
are neglected because of similar energies to the electrons and cannot gain much from 
Compton scattering (Done 2010, Done et al. 2011). We compute the Compton spectrum 
by Monte Carlo simulation. The accretion flow consists of a cold disk in the middle and 
hot corona above/below. At each radius, the cold disk emits black body radiation with the 
temperature T e g as shown in the model equation dH|). The soft photons from the local disk 
transverse through the hot corona. Since the hot corona is optically thin, some of the soft 
photons pass through it without scattering, some of them are scattered by the electrons 
in the corona. We calculate the energy spectra of photons that emerge from the corona. 
The method of Monte Carlo simulation used in our paper is essentially the same as that 
described by Pozdnyakov et al.(1977). In the code, we introduce the weight uj to efficiently 
calculated the effects of multiple scattering. Firstly, we set initial weight = 1 for a given 
soft photon from the cold disk. Then, we calculated the escape probability Pq of passing 
through the corona slab. The value of ujqPq is the transmitted portion and is recorded to 
calculate the penetrate spectrum. The remaining weight u>i = Wo(l — Po) is the portion that 
undergoes more than one scattering. Suppose that P n is the escape probability after the 
nth scattering, then u Q P n is the transmitted portion of photons after the nth scattering. 
The remaining portion u n (l — P n ) undergoes the (n + l)th scattering. This calculation is 
continued until the weight u n becomes sufficiently small. In our calculation, the number of 
processes calculated is 100000 and u n < 10~ 4 . Finally, we obtain the transmitted spectrum 
and the soft spectrum. 



13 



3. NUMERIC RESULTS 

To study AGNs, we fix the black hole mass to be Mbh = 1O 8 M all through our 
calculations and Eddington accretion rate is M Edd = 1.39 x 10 26 gs -1 . The viscous 
parameter a in the corona is fixed to be 0.3. 

3.1. The Properties of Corona under Viscous Heating 

First, we set A = 0.0, meaning / = 0, which is the case that the corona is heated by 
only viscous friction. We calculate the corona structures at R = 1000 R s for rh = 0.02 
and rh = 0.2 as shown in the FigJIJ here rh = M/M^dd which is the total accretion rate. 
Comparing the results in two panels, we find that the vertical profiles of the coronal 
quantities are similar for different accretion rates. From the lower boundary upwards, 
the pressure (P) and the vertical mass flow (pv z ) decrease dramatically in a thin layer. 
The downward heat flux (—F c ) increases in this layer until a maximum is reached and 
then decreases along z and can be neglected on the upper boundary. The temperatures 
of electrons (T e ) and ions (7$), starting at a coupled value of 10 6,5 K at coronal lower 
layer, increase with the height and decouple at about z/R ~ 0.1. Upon this height, 7] still 
increases up to 7] ~ 0.3T v i r and T e keeps almost the same value (~ 0.05T V i r ) throughout 
the corona, where T vir = GMm P n/ kR = 3.37 x 10 12 (i?/i? s ) _1 K is the virial temperature. 
These typical features in disk corona, which are also shown in earlier works (e.g. Meyer 
et al. 2000b), indicate that the corona above a thin disk undergoes very steep changes in 
temperature and density and cannot be simply averaged in the vertical direction. 

We also calculate the evaporation rate rh e va P (= M cvap / M^dd) along distances, which 
represents the mass flowing rate in the corona. The results are shown in Figj2] Similar to 
the previous works (e.g. Liu at al. 2002a), the coronal accretion rate increases toward the 



0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

z/R z/R 

Fig. 1. — Vertical structure of the corona around a black hole at a distance R = 1000 i?s for 
accretion rate rh = 0.02(left panel) and m = 0.2 (right panel). Here, z is the height above 
the equatorial plane of the accretion flow; T e and 7] are the electron temperature and ion 
temperature, respectively, while T v ; r is the virial temperature used for scaling. F c , pv z , P, 
and v z /v s are the heat flux, vertical mass flow rate per unit area, total pressure, and vertical 
velocity scaled by the sound speed, respectively. The quantities at the lower boundary are 
marked with subscript 0. 
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central black hole, reaches a maximum at a few hundred Schwarzschild radius and then 
drops very quickly in the inner region because of strong inverse Compton scattering, i.e. a 
large fraction of the corona gas condensing onto the disk near the black hole. Comparing 
the three evaporation curves for accretion rate, m=0.02, 0.2 and 0.4, we find that the 
maximum evaporation rate decreases with increasing accretion rate and the corona gas 
begins condensing at outer regions for higher accretion rate. The temperature in electrons, 
as shown in Fig.l, is lower at higher accretion rates. This indicates that the corona becomes 
very weak at high accretion rates, like in HLAGNs. 

-1.5 
-2.0 

a 
o 

-3.0 
-3.5 

2.5 3.0 3.5 4.0 4.5 

lg(R/Rs) 




Fig. 2. — The relation between the corona accretion rate m ev ap and the radius for different 
accretion rate with / = 0.0. 



We can roughly estimate the radiation from the disk and corona under the viscous 
heating. As shown in Figj2], the maximum corona accretion rate is m e vap,max = 4 x 10~ 3 at 
~ 1250 R s for rh = 0.2. Supposing that the corona keeps accreting at this accretion rate, 
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the radiation strength of corona relative to the disk is 

A: _ ?]M C C ^ ^evap.max _ q Q2 (^9) 

L d r]M d c 2 ~ rh- rh evap 

.max 

Here, we have supposed that radiative efficiency r] in corona is equal to that in disk. From 
Eq. (]19p . we can see that the corona is very weak and the radiation is dominated by the disk. 
In fact, the luminosity ratio between the corona and disk is much smaller than above value 
since more and more coronal gas condenses into the disk at smaller distance (R < 1250i? s ). 
In order to produce strong X-ray in corona, we consider that some fraction of gravitational 
energy is transferred to the corona, providing additional heating mechanism against the 
strong Compton cooling. 



3.2. The Properties of Corona with Additional Heating 

3.2.1. Vertical structure 

Assuming a fraction / (=0.1, 0.2, 0.3 and 0.4) of accretion energy is added to the 
corona, we perform numerical calculations for accretion rates of rh =0.1, 0.2 and 0.5. 

We start our calculations from the distance of 100i? s to the ISCO region (R = 3.1i? s 
for convenience). It has been shown in section [37TI that the maximum corona accretion 
rate can only reach ~ 4 x 10~ 3 MEdd at ~ 1000i? s for rh = 0.2. In the inner region, 
strong Compton cooling leads to condensation of corona. Whether the corona can exist at 
distances inside 100R S depends on how much accretion energy is added to the corona. From 
our test calculations, we find that for / = 0.2, a corona can still survive in the region from 
100i?s ~ 3.1.Rs- The coronal vertical structure is shown in Fig]3j The profile is similar to 
the case for / = at large distances, though the absolute values are different. 

In order to study the corona properties with different energy input, we calculate the 
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Fig. 4. — The radial distribution of corona accretion rate for rh = 0.2 (left) and m = 0.5 
(right) with different /. 
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coronal accretion rate for / = 0.1, 0.2, 0.3, 0.4, where the total mass supply rates are 
assumed to be m = 0.2, 0.5 respectively. The results are plotted in FigJUand FigJSJ One 
can see from FigJUthat the accretion rate in the corona decreases with decreasing distance, 
as a consequence of condensation. While the electron temperature and scattering optical 
depth are nearly the same at distances up to 100i? s , as shown in Fig. El Besides, Fig. HJalso 
shows that the larger fraction / is, the larger becomes the gas flow in the corona. This is 
because more energy goes directly to heat the electrons in corona to a higher temperature, 
leading to more conductive flux and hence more evaporation or less condensation. Thus, 
more gas is kept in and accreted through the corona. This is confirmed by the higher 
electron temperature and larger scattering optical depth for larger /, as shown in Fig. |5j 
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Fig. 5. — The radial distribution of scattering optical depth (r) and electron temperature 
(T e ) for rh = 0.2 with different /. 

Comparing with the two panels in FigHJ we understand that the coronal accretion 
rates are insensitive to the mass supply rate, as long as the fraction of accretion energy 
released in the corona is the same. Calculations show that the electron temperature and 
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scattering optical depth of the corona do not vary with the accretion rate either. This can 
be understood as a consequence of energy balance in the corona. For given /, when the 
mass supply rate (m) increases, the heating rate added to the corona also increases (oc ffn). 
On the other hand, the cooling rate by inverse Compton scattering also increases with the 
increasing density of soft photon (oc (1 — f)rh, if (1 — f)rh » m cvap ). If the additional 
heat and Compton cooling are the dominant energy terms for the coronal electrons, there 
should be no variation in both the evaporation rate and Compton ^-parameter. Therefore, 
the electron temperature does not vary either. 

3.2.2. Spectra of the disk and corona 

Deriving the structure of our model at each radius, we get radius- dependent electron 
temperature and electron scattering optical depth, which are the input parameters for 
bremsstrahlung and Compton spectrum. The bremsstrahlung radiation is proportional 
to nl but the inverse Compton scattering is proportional to n c . So from the vertical 
structure, we can see that the bremsstrahlung dominates in the lower layer of the corona 
and the upper layer is dominated by the Compton radiation. Since the temperature and 
number density in the upper corona still depend on the vertical height, we need to choose 
appropriate temperature and density at every given distance for calculating the Compton 
spectrum with Monte Carlo simulation. Here we constrain them by the luminosity, i.e., 
when the luminosity derived from the Monte Carlo simulation equals to the one derived 
from the coronal structure calculation, the chosen temperature and density are thought to 
be right and we take this spectrum as the true spectrum from the corona. 

FigE] shows the spectra for different f,f = 0.1, 0.2, 0.3, 0.4, and different accretion 
rates, rh = 0.2, 0.5. The spectrum is Compton dominated in the high accretion rate. The 
larger / is, the stronger becomes both Compton and bremsstrahlung radiation, since the 
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Fig. 6. — The bremsstrahlung and Compton spectrum for rh = 0.2 (left) and rh = 0.5 (right). 
The curves with the bump in about 10 19 Hz respond to the bremsstrahlung emission, and 
the ones with bump in the UV range are contributed by the blackboday radiation in the 
disk and inverse Compton scattering in the corona. The two vertical dotted lines indicate 
the photon energies of 2 keV and 10 keV, respectively. 
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Fig. 7. — The total spectrum contributed from multi-blackbody from disk, Bremsstrahlung 
and inverse Compton scattering for / = 0.2 (left panel), / = 0.4 (right panel). 

coronal density and temperature increase with increasing heating to the corona; While the 
optical-UV radiation, which is from the disk, slightly decreases with increasing /. The 
spectrum of hard X-ray in the range of 2 — 10 keV is marked by the two vertical dashed 
lines. The hard X-ray spectrum is harder for higher / and L^, \/ L 2 -\q kcv decreases with /. 
For / = 0.1, the spectrum seems so steep that the corona radiation is still very weak. The 
photon index for the hard X-ray is in the range of 2.2 — 2.7 for 0.2 < / < 0.4, which roughly 
fits the spectrum of HLAGNs. Table [1] lists photon index for 2 — 10 keV and the ratio of 
-^boi/-^2-io kcv- F° r different accretion rates, the spectrum index is nearly the same with 
the same /, which can also be seen from the total spectrum shown in Fig. [3 Comparing 
with the two panels in the FigJTJ we understand that the spectrum varies with the fraction 
of accretion energy liberating in the corona, while it does not sensitively depend on rh, as 
long as / is independent on the accretion rate. This can be understood from our structure 
calculating results, namely, for given / the electron temperature in corona T e , and the 
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scattering optical depth r, do not change as the accretion rate, though the disk radiation 
increases and hence Compton radiation also increases with increasing accretion rate. 

The photon index and the luminosity ratio, as listed in Table 1, indicate that a 
fraction of accretion energy is necessary in order to explain observations in HLAGNs. For 
0-2 < / < 0.4, the predicted photon index is around 2.2 — 2.7, and the bolometric correction 
from hard X-rays Lboi/^2-10 koV is in the range of 28 — 228. These are roughly consistent 
with observations (Vasudevan & Fabian 2007, 2009; Shemmer et al. 2006 and references 
therein; Zhou & Zhao 2010). In Fig. E]we show how the photon index in X-rays varies 
with bolometric luminosity. The figure reveals that the bolometric correction is larger for 
a steeper X-ray spectrum, which seems to be a common feature for objects with different 
Eddington ratios. This is a consequence of decrease of energy fraction released in the 
corona. We now raise the questions as, what is the underlying physics driving the accretion 
energy to release in the corona? Does the energy fraction / depend on the accretion rate or 
the black hole mass? We discuss this in the following section. 

Table 1: The photon index and L bo i/L 2 _ 10 kcV for different accretion rate with different /. 
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Fig. 8. — The relation between photon index (T ) and bolometric correction (Lboi/^-io kev)- 
It shows that T increases with Lboi/-^2-io kcV) which is a consequence of decreasing energy- 
fraction released in the corona. 



-24- 



4. DISCUSSION 

4.1. Heating Mechanism in the Corona — Magnetic Heating? 

It has been a long history in studying magnetic field in the accretion flow, either with 
magnetohydrodynamics(MHD) simulations (e.g. Balbus & Hawley 1998; Matsumoto 1999; 
Stone & Pringle 2001; Machida et al. 2001; Machida & Matsumoto 2003; Kato et al. 2004; 
Kawanaka et al. 2008 and references therein; Ohsuga, et al. 2009, 2011; Penna et al. 2010) 
or numerical simulation to interpret observations (Liu, et al. 2002b, 2003; Cao 2009). 

There are good reasons to believe that the magnetic field plays an important role in 
transporting the accretion energy from the disk to the corona (Haardt & Marashi 1991). 
Liu et al. (2002b; 2003) investigated the Parker instability in the accretion flow, following 
the model for the solar corona (e.g. Yokoyama & Shibata 2001) and found that the 
accretion energy can be transported to the corona and heats the electrons in the corona. 
In the frame of disk corona with mass evaporation, Qian, Liu & Wu (2007) analyzed the 
coronal structure and evaporation features with the magnetic field and found that the 
maximum evaporation rate stays more or less constant (~ 3% Mgdd ) when the strength of 
the magnetic field changes, but that the radius corresponding to the maximum evaporation 
rate decreases with increasing magnetic field. Qiao & Liu (2009) investigated the influence 
of viscosity parameter a in the corona and suggested that the coronal radiation is largely 
enhanced when viscosity increases. 

However, it is not very clear how magnetic field connects with viscosity in the 
evaporation-fed corona, though some investigations of accretion flows show that the 
magnetic field and viscosity are associated (e.g. Balbus & Hawley 1991; Matsumoto & 
Tajima 1995; Abramowicz et al. 1996; Balbus 2003). Here we discuss whether the predicted 
spectrum can be consistent with observations if it is the magnetic field that brings a fraction 
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of gravitational energy to the corona. 

Let us assume that the magnetic field is generated in the disk by dynamo action. 
Then the magnetic loops can emerge out of the disk by Parker buoyant instability. In this 
process, the accretion energy stored in the magnetic field is transported into the corona and 
thereby liberated through magnetic reconnection. According to the equipartition theorem 
between the magnetic pressure and gas pressure, we set Pb = f^- = i-fg,d, where Pb and 
P gi d are magnetic pressure and gas pressure in the disk. Then the magnetic energy flux is 



If this part of magnetic energy flux is the additional heating source for the corona, the 
ratio of this part of energy to total gravitational energy is 

f = TT-= i-a P ,A) l P7lF t0 u (21) 

-Ttot P 

where F tot is the total accretion flux through accretion process, F tot = 3 ^^ 7 [l — (3P S / P) 1 / 2 ] 
8.56 x 10 26 mm~ 1 r~ 3 $, where $ = 1 - (I) 5 . 

In the thin disk, the pressure of the disk is 

P = P g ,d + Pb + P,d = (1 + I) ^ + ^£ (22) 

p fim p 3 

where \i = 0.62 is the averaged molecular weight. Since some fraction (/) of gravitational 
energy is carried into the corona by magnetic activity, then the energy equation in the disk 
is 

= (1 - /) x F tot = (1 - f)-^~ [1 " (3P S /P) 1/2 ] , (23) 

where r = (n es + K p d T d " 3 ' 5 )E, n cs = 0.4 cm 2 g _1 , k = 6.4 x 10 22 cm 2 g" 1 and S = 2p d H d 
is the surface density of the disk. And the angular momentum equation in the disk is 

M 

-$ = i/E, (24) 

07T 
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where v = ^ac s H d , the viscous velocity c s = QH d . 

If the gas pressure in the disk is dominant over radiation pressure, P gi d 3> P r ,d, and 
Thomson scattering is much higher than free-free absorption, ^> erg, we have r ~ K es S) 
which is similar to the case b) in Shakura & Sunyeav (1973) except that r = R/R s is taken 
here instead of r = P/3P S . We can derive the disk quantities. The density in the disk is 

p A = 21.43(am)~^ [m$]* r~i(l + 77)^(1 - f)~^, (25) 

and the gas pressure is 

P g , d = 1.91 x 10 18 (am)"^ [m$]*H&(l + 4)~s(l - /T^- (26) 

With this disk pressure and density, the energy fraction can be obtained from Eq. (l2ip . 



/ = 1.88a- 1 (l + / 9)-i. (27) 

One can see that / is independent on the accretion rate if the disk is gas pressure dominated. 
This implies that a constant fraction of disk accretion energy is transported to the corona 
at different accretion rates, if equipartition coefficient (3 does not vary with accretion rates. 
If this is the case for the additional heating to the corona, the model spectrum will not vary 
with the accretion rate, as shown in FigJTl 

In fact, the accretion disk around a supermassive black hole is often dominated by 
radiation pressure. In the case of P gj d <§C P r .d and ctt S> crs, the density and pressure in the 
disk are expressed as, 

p d = 1.14 x 10~ 6 (am)' 1 [m$]" 2 ri(l - /)" 3 , (28) 
P gid = 6.21 x 10 9 (am)"t [m$]" 2 ri(l - f)~^, (29) 
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which are similar to the solution of case a) in Shakura & Sunyeav (1973). Then we can 
determine / from Eq. f[2"Tj) as 

/(l - /)¥ = 1.52 x 10~ 9 a-^m-i(m$)- 3 ri/3-t. (30) 

Apparently, in a radiation pressure dominated disk, / depends on the accretion rate. 
The function /(l — /)x in Eq. (130]) reaches a maximum of 0.095 at / = Jr, which requires 
1.52 x 10~ 9 a - ^m~l(m$) -3 rT§/3 - 5 < 0.095. This indicates that the accretion rate must be 
larger than a certain value in order to get a solution for / from Eq. ( 1301) . If this condition 
1.52 x 10~ 9 a~^m _ l(m$) _3 rii/3 _ f < 0.095 is fulfilled, there are two sets of solutions for / 
( We saw this bimodal nature in Liu et al. (2002b)). As shown in Figj9j one set of solutions 
is in the range of < / < Jg, where / decreases with increasing accretion rate rh; While the 
other one is in the range of Jr < / < 1, where / increases with increasing rh. The relation 
between the energy fraction and the accretion rate, i.e. /(rh), depends on the viscosity 
and magnetic field. In the left panel of FigJHl we show the effect of viscosity by assuming 
= 1 under the theory of equipartition of gas pressure and magnetic pressure. Actually, 
the local and globe MHD simulation in radiation-dominated accretion flow reveals that 
magnetic energy density may exceed gas energy density but kept below radiation energy 
density, which means that (3 can be less than 1.0 (e.g. Turner et al. 2003; Ohsuga et al. 
2009, 2011). The effect of magnetic fields is plotted in the right panel of Figj9j It shows 
that the influence of magnetic field is significant. The stronger the magnetic field, the larger 
rh corresponding to / = 8/35. For a given accretion rate, for example, rh = 0.4, / is larger 
with stronger magnetic field (i.e. smaller /3) in the lower part of the curve (/ < 8/35). 

Observed spectral features provide a clue to possible range of /. If / < Jg, the 
fraction of accretion energy transferred to the corona decreases with accretion rate. As a 
consequence, the hard X-ray spectrum becomes softer at high accretion rates, as shown in 
FigJTOl However, if a large fraction of accretion energy (/ > Jr) is transferred to the corona, 
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Fig. 9. — The fraction of accretion energy released in the corona as a function of accretion 
rate. For a given accretion rate, there are two values for /, one is / > ^, the other is / < Jr. 
The horizontal line indicates / = 8/35, above which / increases with increasing accretion 
rate and below which / decreases with rh. The effect of viscosity is shown in the left panel 
by assuming (3 = 1 and the effect of magnetic field is shown in the right panel by assuming 
a = 0.01. It can be seen that for small viscosity and/or strong magnetic field, the accretion 
rate corresponding to a certain / is large. 
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/ increases with accretion rate, consequently, the spectrum is harder at higher accretion 
rates. Comparing the model prediction and observations, a small fraction of accretion 
energy liberation in the corona is preferred, that is, < J|. In FigJlOj we plot the spectrum 
for different accretion rate. The dashed lines represent a spectrum for a low accretion rate, 
and the dotted line is for a relatively higher accretion rate, and the solid line is for a much 
higher accretion rate, where the fraction of energy released in the corona corresponds to 
0.3, 0.2, 0.1 respectively, determined by Eq.f l30|) assuming disk viscous parameter a = 0.01, 
= 1 for a m = 10 8 black hole at the most efficient radiation region r = 49/12. We 
show that the spectrum becomes steeper/softer at high accretion rates. Consequently, the 
bolometric correction from the 2-10 keV increases with increasing accretion rate. 

We point out that, in the high accretion rate, the radiation pressure dominated disk 
is thermally and viscously unstable, which can lead to the S-shaped "limit-cycle" (e.g. 
Shakura & Sunyaev 1976; Abramowicz et al. 1988; Honma, Matsumoto & Kato 1991; 
Szuszkiewicz & Miller 1998; see however, Hirose et al.(2009) who claimed no instability 
by local radiation- hydro dynamical simulations). However, if a fraction of accretion energy 
is liberated in the corona or the outflow the disk can be stable (Nakamura & Osaki 1993; 
Svensson & Zdziarski 1994). The magnetic heating to the corona discussed here is helpful 
to stabilize the disk, though / is not large enough. 

4.2. Effect of viscous heating directly to electrons 

In our corona, we assume that all of the viscously dissipated energy heats ions only. 
However, direct heating to the electrons can also increase the corona emission. According 
to the detailed modeling to Sgr A* (Yuan et al. 2003), a significant fraction of the viscously 
dissipated energy should heat electrons directly. Such a result was later supported by the 
numerical simulation done by Sharma et al. (2007). Nevertheless, this effect is negligible if 
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lg(i//Hz) 



Fig. 10. — The spectrum for different accretion rate. The dashed lines represent a spectrum 
for a low accretion rate, the dotted line indicates a higher accretion rate, and the solid line is 
for much higher accretion rate. The fraction of energy released in the corona for theses three 
curves corresponds to 0.3, 0.2 and 0.1, respectively. The corresponding m is determined by 
Eq.( )30|) by assuming disk viscous parameter a = 0.01, /3 = 1 for m = 10 8 black hole at the 
most efficient radiation region r = 49/12. 
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much larger additional heating to the corona is assumed. When a strong disk exists under a 
corona as in luminous AGN, a small fraction of disk accretion energy transported to coronal 
electrons can dominate the coronal heating. For instance, assuming m = 0.2 and / = 0.1, 
the additional heating to the corona is 0.02 (in unit of Eddington luminosity); While the 
accretion rate in the corona, as shown in Fig. 4, is 10 -4 to 10~ 3 in the major radiation 
region, which produces viscous heat of 10~ 4 to 10~ 3 (in unit of Eddington luminosity). Even 
if all of this viscous heating directly goes to electrons, it is much smaller than the additional 
heating 0.02. The ratio is in the range of 0.5% to 5% . This ratio is nearly the same for a 
higher / since the accretion rate in the corona increases with increasing /. For an accretion 
rate higher than 0.2, the ratio is smaller, which can be derived from the accretion rate of 
corona shown in the right panel of Fig. 4. Therefore, the neglect of viscous heating directly 
to electrons is a reasonable approximation. 

4.3. Disk Wind at High Accretion Rate 

When the accretion is greater than O.lMEdd, the line driven wind and radiation 
pressure driven wind will be important, which results in the accretion rate decreasing with 
the decrease of radius as M oc R a (0.0 < a < 1.0) (see Proga et al. 2000; Ohsuga et al. 
2005). It is clear that with wind loss the soft-photon field becomes weaker than the case 
without winds. Therefore, the Compton cooling is weaker, leading to less condensation. 
When the system reaches equilibrium, the spectrum can be harder. This is the case without 
additional energy input. In this work, an additional energy is added into the corona and 
we find that the larger the / is, the harder the spectrum is. Since < a < 1, the energy 
generating rate oc GA, ^ M oc is still dominated by the innermost region. Thus, the wind 
loss causes a decrease of disk luminosity and hence a harder spectrum. However, if wind 
mass loss is very significant, there could be no longer optically thick disk and the effects are 
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large. 



4.4. Irradiation 

A corona lying above a cool disk can illuminate the disk, which increases the energy 
density of soft photons, leading to stronger coronal radiation. This effect can be significant 
in the case of a strong corona above a truncated disk. However, for a relatively weak corona 
above a strong disk, as the case of HLAGNs, the effect of irradiation is unimportant. We 
estimate this effect as follows. 

Photons emitted from corona are roughly isotropic, about half of which go down to a 
slab-like disk, irradiating the disk. The irradiating flux is, 



F in = \f c {1 -a)= l -fF tot {l - a), (31) 

where a is albedo which is usually taken as 0.15 (e.g. Zdziarski et al. 1999). This flux is 
scattered/absorbed in the disk and eventually re-emitted, contributing to the soft photon 
field. The ratio of this illuminating flux to the disk radiation, = -F tot (l — /), is 



F irr _ |/(1 - a) 



(32) 



F d (1 - /) 

This ratio reaches 1 when / > 0.7, indicating that the irradiation becomes dominant only 
when a large fraction of accretion energy is released in the corona. Therefore, we don't 
include the irradiation in this work in order to make the numerical calculations simple. 



4.5. Comparison with the MHD Results 

Hirose, Krolik & Stone (2006) performed the shearing box numerical simulation 
with radiative transport, focusing on the vertical structure of the standard thin disk 
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around 300i? g (R g = ■ They found that there is actually very little dissipation in 
the coronal region as the corona is magnetically dominated. Consequently, the corona 
is too weak to explain the hard X-ray emission observed in AGNs. The disk corona 
evaporation/condensation model without additional heating gives similar results as MHD 
simulation (for details see Meyer-Hofmeister, Liu, & Meyer 2012). When a corona is heated 
only by viscous dissipation within the corona, efficient inverse Compton scattering of the 
disk photons leads to over-cooling of the corona. As a consequence, part of the coronal gas 
condenses into the disk, leaving a very weak corona. 

To prevent the corona from collapsing, we assume a fraction of accretion energy is 
transported from the disk and release into corona in this paper, which appears conflict 
with the MHD simulation results. We note that in MHD simulation of Hirose et al. 
(2006), it was presumed that the disk is dominated by gas pressure, which is reasonable 
for studying the vertical structure of black hole X-ray binaries in about 300i? g . However, 
we are investigating the inner region of AGN disk, where the disk is radiation-pressure 
dominant. Recently, Blaes et al.(2011) simulated such kind of accretion flow and found that 
the thermodynamics of a radiation-pressure dominant disk differs significantly from that 
envisaged in standard static models of accretion disks. In radiation-dominant plasma, the 
buoyant motions become significant for the overall energetics of the plasma. The photons' 
outward advection becomes comparable to radiative diffusion, and the associated vertical 
expansion work balances vertical heat transport and dissipation. At the present stage, it is 
not clear how much discrepancy exists between MHD simulation and our results. 



4.6. Luminous Hot Accretion Flow as an Alternative Model 

Based on the thermal stability analysis presented in Yuan (2003), luminous hot 
accretion flow (LHAFs, Yuan 2001) could also have a two-phase structure. Yuan & 
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Zdziarski (2004) calculated the emitted luminosity from such kind of two-phase accretion 
flow and found that it can produce X-ray luminosity as high as 10% Eddington luminosity 
This is an alternative model for AGNs with intermediate spectrum and Eddington ratio. In 
our work, we aim at more general case for HLAGNs with luminosity up to Eddington value, 
and calculate the broad waveband spectra from both the disk and corona. 

5. Conclusion 

We applied the disk evaporation model to the high-luminosity AGNs. To explain the 
typical spectra of HLAGNs, we conclude that there should be additional heating to the 
corona to prevent it from over-cooling by strong inverse Compton scattering. 

We assume that a fraction of gravitational energy (/) is liberated in the corona and 
calculate the corona structure. Then the spectrum from the disk and corona is calculated 
by Monte Carlo simulation. We find that the hard X-ray is dominated by Compton cooling 
and the photon index for hard X-ray in 2 — 10 keV is 2.2 < T < 2.7. For a given accretion 
rate, T and £ft //-^2-iokcv decrease with increase of /. We discuss a possible mechanism for 
the presumed heating to the corona, that is, the magnetic heating. With equipartition of 
magnetic energy to the gas energy in the disk, we derive the fraction of accretion energy 
stored in magnetic field. Assuming this energy is released in the corona, the model predicts 
that, for / < Jjr, the hard X-ray becomes softer at higher accretion rate and L^i /-^2-iokev 
increases, which is roughly consistent with the observational results. 
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